Coupling between radial and torsional oscillations 
in a magnetized plasma and possible stellar applications 
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ABSTRACT 

The weakly non-linear interaction of slow magnetoacoustic and torsional Alfven waves prop- 
agating along an unperturbed poloidal magnetic field in a stellar interior is studied in the high 
plasma /? limit. It is shown that slow magnetoacoustic waves parametrically drive torsional modes 
with the half frequency and wave number. Therefore global slow magnetoacoustic waves which 
have a radial velocity polarization, may amplify the torsional oscillations. Possible applications 
of this mechanism to the Sun and binary stars are briefly discussed. 

Subject headings: Starsimagnetic fields - Stars:oscillations - Magnetohydrodynamics (MHD) - Waves 



Introduction 



Al though Walen's original suggestion (jWalen 
1949) that torsional azimuthal Alfven oscillations 
may exist in stellar interiors, the energy source 
sustaining these oscill ations is not yet known 
( Rosner fc Weisd Il992ft . Without an energy 
source, the oscillatio ns are believed to be damped 
by p hase mixing (jCharbonneau fc MacGregor 
19931 ) or transport of magnetic flux towards the 
surface. Several different mechanisms have been 
suggested fr om time to time to s upport torsional 
oscillations (lLavzer et aL 19791 : Mclntvre 1994 ; 
Lanza et al.lll998t lGoughll2002l ). however none of 
them may account for the energy source. The 
main problem is that torsional Alfven waves are 
pure electromagnetic oscillations in magnetohy- 
drodynamics and one needs to have an azimuthal 
force which may support them against damping. 
The absence of such a force suppressed the inter- 



est in these oscillations as a possible mechanism 
for stellar activity cycles. However no one paid at- 
tention to the fact that it is not necessary to have 
a direct driving force for the amplification of the 
oscillations. The oscillation can also be amplified 
through the periodical variation of the system pa- 
rameters and the process is called parametric res- 
onance. The well-known example of such system 
is the mathematical pendulum with periodically 
varying length. When the frequency of the param- 
eter variation is twice the frequency of the system 
oscillation, then the oscillations will exponentially 
grow in time. The phase speed of torsional Alfven 
waves depends on the magnetic field and medium 
density. So if some external force (or oscillation) 
causes a periodical variation of either quantities 
(or both together) then the Alfven waves will be 
ampli fied. The pioneering work of Zaqarashvilil 
(|2001l ) suggested that the periodical density vari- 
ation caused by sound waves in the high plasma (3 
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limit may amplify the s h ear A lfven waves. Then, 
Zaqarashvili fc Roberts! (|2002h showed the reso- 
nant energy transfer from fast magneto-acoustic 
waves propagating across the unperturbed mag- 
netic field to Alfven waves propagating along the 
field. The energy transfer is relative to the veloc- 
ity component of resonant Alfven waves normal to 
both magnetic field and direction of propagation 
of the magnetoacoustic waves. In both papers the 
temporal behaviour of Alfven waves was found 
to be governed by Mathieu equation. Thus the 
frequency of resonant Alfven waves was the half 
frequency of the driving compressible oscillations. 
The natural astrophysical suggestion was that 
stellar pulsation in the fundamental mode may 
ampl ify torsional Alfven wave s in a seed magnetic 
field IjZaqarashvili et al l2002h . The process can be 
of importance in stellar chromospheres and coro- 
nas as the Alfven waves propagate upwards carry- 
ing energy and momentum and so may contribute 
to chromospheric and coronal heating. Also the 
waves m ay propagate up to the stel l ar win d and 
recently IZaqarashvili and Belvedere] ( 2005 ) sug- 
gested them as a source for long period Alfven 
oscillations (with periods of few hours) observed 
by ULYSSES in the solar wind. On another hand, 
the fundamental frequency of stellar pulsation is 
too high for resonant torsional Alfven waves to 
give rise to large-scale torsional oscillations in a 
stellar interior. Therefore one needs a longer pe- 
riod pulsation in order to suggest them as an 
energy source for torsional oscillations. 

Here we suggest that the slow magnetoacoustic 
waves propagating along the magnetic field in the 
high plasma [3 limit are nonlinearly coupled to the 
torsional Alfven waves propagating in the same 
direction. The velocity polarization of slow mag- 
netoacoustic waves propagating along a straight 
uniform magnetic field is along the radial direc- 
tion in the case of a cylindrical coordinate sys- 
tem (the waves have a velocity component along 
the magnetic field as well). Therefore fixing the 
boundaries along the axis of the cylindrical sys- 
tem originates a standing pattern of waves lead- 
ing to the pulsation of the magnetic cylinder with 
a spatial scale which corresponds to the first or 
second eigenmode. In the case of stars with a 
poloidal magnetic field in the interior, the funda- 
mental mode of slow magnetoacoustic waves has a 
very low frequency due to the small phase speed 



(note, that the slow magnetoacoustic waves prop- 
agate with the Alfven speed in high (3 limit). So it 
will lead to long period stellar pulsation. Since the 
velocity of the slow magnetoacoustic waves has the 
radial direction, any force (gravity, radiation pres- 
sure etc) may support this oscillation. If the oscil- 
lation transfers energy to torsional Alfven waves, 
then they will lead to the set up of torsional os- 
cillations, as their wavelength will comparable to 
the stellar radius. So the torsional oscillations will 
gain energy from the radial pulsation and thus can 
be sustained against damping. Here we consider a 
simple case of wave coupling in an uniform mag- 
netic cylinder, which can be accounted as the sim- 
plest model of a non-rotating star with a poloidal 
magnetic field. 

2. Equations and mathematical develop- 
ments 

We consider a non-turbulent medium with zero 
viscosity, infinite conductivity and negligible dis- 
placement current, which can be described by the 
magnetohydrodynamic (MHD) equations: 



dB 

~at 



Vx(VxB) 



(1) 



P^+p(V-V)V = -VP+i(VxB)xB, (2) 
at fx 



dp 
dt 



V(pV) = 0, 



(3) 



where p is the medium density, P is the pressure, 
V is the velocity, B is the magnetic field and // is 
the magnetic permeability. To this system we add 
the adiabatic equation of state 



P = Po 



(4) 



where po and po are the unperturbed pressure and 
density, while 7 is the ratio of specific heats. We 
adopt a cylindrical coordinate system (R, 4>, Z) 
and for simplicity consider only the axisymmetric 
problem, so everywhere d/d(j) = is assumed. We 
suppose that an uniform magnetic field is directed 
along the Z axis so B — (0, 0, Bo). 

The medium is supposed to be uniform, there- 
fore the equilibrium is satisfied automatically. It 
must be mentioned, however, that real astrophys- 
ical situation is much more complicated. In the 
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stellar interior the gravity and the pressure gradi- 
ent necessarily lead to the non uniform spatial dis- 
tribution of unperturbed physical quantities (den- 
sity, temperature and magnetic field). But in or- 
der to understand the physics of the wave coupling 
we choose the simplest uniform model (Fig.l), 
which considers an uniform magnetic cylinder with 
an uniform density. Thus, gravity is neglected in 
the momentum equation (note that almost all dy- 
namo models also ignore the gravity). Once the 
physical grounds of the coupling are clear, future 
models may include a more realistic spatial distri- 
bution of the physical variables. 

For the study of the coupling between different 
oscillations we use the perturbative method i.e. all 
variables are perturbed around their equilibrium 
values and their consequent evolution is studied. 

We first briefly describe the linear evolution of 
the perturbations in form of linear waves and then 
study the weakly non-linear interaction between 
the waves, which leads to the coupling of our in- 
terest. 

3. The linear approximation 

All physical variables are presented as a sum 
of unperturbed and perturbed parts ^ = Co + Cl- 
in the linear approximation only the first order 
terms are retained. This allows us to consider the 
evolution of different waves separately. The sys- 
tem (l)-(4) includes Alfven and magnetoacoustic 
waves. 

3.1. Torsional Alfven waves 

In cylindrical coordinates the pure Alfven waves 
are the torsional waves, which have azimuthal ve- 
locity and magnetic field components and repre- 
sent the purely electromagnetic part of oscilla- 
tions. They are described by equations 



914 



(5) 



(6) 



at " u dz ' 
Po lk ~ ~~dz' 

where V&, are the velocity and magnetic field 
perturbations. As we already mentioned, the un- 
perturbed density is considered to be homoge- 
neous. 

The restoring force for the torsional waves is 
the tension of magnetic field lines only. The waves 



propagate along the magnetic field lines with the 
Alfven speed 

CA = ~7==- ( 7 ) 



Consequently the dispersion relation of the 
waves is 



(8) 



where lu a is the frequency and k A is the wave num- 
ber. 

The waves are strictly transversal and they do 
not cause a density perturbation in the linear 
regime. 

3.2. Slow magnetoacoustic waves 

We consider the adiabatic process, then the re- 
maining linearized equations are 



Po 



db R _ 8V R 

db z = Bp d{RV R ) 

dt R dR ' 

dV R dpi B f db R _ dbz_ 

dt ~ dR [i \dZ dR 

dVz_ dpi 
~~ dZ' 



po 



dpi 
dt 



-Po 



dt 

1 d{RV R ) dV z 
R dR dZ 

Pi = c 2 s pi, 



(9) 

(10) 

(11) 

(12) 

(13) 
(14) 



where V R , Vz, b R , bz, pi and pi are the perturba- 
tions of velocity, magnetic field, density and pres- 
sure respectively; c s = y jpp Jpo is the adiabatic 
sound speed. 

After Fourier analysis with ~ exp (— iujt + ik z Z) 
Eqs. (9)- (14) lead to the equation 



d 2 V R 1 dV R 



dR 2 R dR 



R? 



V R = 0, (15) 



where 



u 2 (cl+c A )-c A c 2 s kl ■ 



(16) 



Eq. (15) is the Bessel equation when m 2 > 0, the 
modified Bessel equation when m 2 < and Eu- 
ler eq uation when m = (jAbramowitz fc Stegun 
19641 ). Here m plays the role of the radial wave 
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number. For m^O the expression (16) leads to the 
dispersion relation of obliquely propagating mag- 
netoacoustic waves. For the parallel propagation 
(i.e. m = 0), the waves propagate either with 
sound or with Alfven speeds; for the high (3 plasma 
(c s va) the fast magnetoacoustic waves propa- 
gate with the sound speed (u = ±c s k z ), while the 
slow magnetoacoustic waves propagate with the 
Alfven speed (lu — ±CAk z ). 

For m 2 > the solution to Eq. (15) can 
be expressed in form of t he Bessel functions 
(jAbramowitz & Stegunlll964l ) 



dependence on R due to Eqs. (9) and (19) 

v R h h R 

Vr = u R —, b R = tin- 
no Kq 



(20) 



while bz, Vz, Pi remain independent on R. The 
wave dispersion relation is the same as for the tor- 
sional waves. The relation between the different 
variables in the slow magnetoacoustic waves prop- 
agating along the magnetic field is as follows: 



Vr = cxJximR) + c 2 Y 1 (mR), 



(17) 



Pi = -2i 



V z = -2% 



p uj 



■UR, 



(21) 



where c\ and C2 are constants. 

For the parallel propagation (m = 0) the solu- 



riR = u R , b z = -2i——ur. 

UJ Ll>Rq 



tion to Eq. (15) can be expressed as (jAbramowitz fc 
19641 ) 

V R = ciR+%, (18) 



where ci and ci are constants. This is the solution 
to slow magnetoacoustic waves. For the fast mag- 
netoacoustic waves Vr = 0. Imposing that the 
velocity is finite on the axis we obtain 



Stegurj t jg easy to snow from Eqs ( 14 ) and pi), that the 
total pressure perturbation at the surface vanishes 
for slow magnetoacoustic waves in the case of high 
[3 plasma (c s ^> ca) due to the dispersion relation 
uj = ±CAk z . Indeed we have for the total pressure 
perturbation: 



Pi 



Bnb 



0«z 



Vr = Cl R. 



(19) 



2i——ur — 

RqUJ \ KJ 



o{(}- 1 ) 



This solution is not valid for unbounded medium 
as it diverges for large R. Thus in cylindrical 
coordinates for unbounded medium slow magne- 
toacoustic waves can not propagate strictly along 
the magnetic field (unlike to the Cartesian co- 
ordinates). However as far as we are concerned 
in lower order harmonics of oscillations with the 
wavelengths comparable to the stellar radius, then 
it is unjustified to consider unbounded medium 
in radial direction. Therefore we should give the 
boundary conditions at the stellar surface. Due 
to the sharp density jump at the surface we may 
set up free boundary conditions, for which the ex- 
pression (19) is valid. Then the velocity has, at 
the surface, the value Vr = c\Rq, where Rq is the 
radius of the star. 

For simplicity, we consider the slow magnetoa- 
coustic waves propagating along an unperturbed 
magnetic field in the medium bounded along ra- 
dial direction. Then the free boundary conditions 
at the surface allow the radial components of ve- 
locity and magnetic field Vr, bR. to have the linear 



therefore the boundary condition is satisfied in the 
zero order approximation with respect to /3 _1 . 

Thus the slow magnetoacoustic waves lead to 
the periodical modulation of the Alfven speed due 
to the variation of density and magnetic field. 
Therefore they may cause a nonlinear influence on 
torsional Alfven waves. In the next section we 
study the phenomenon in the weakly non-linear 
approximation. 

4. The non-linear coupling 

The non-linear influence of slow magnetoacous- 
tic waves on the dynamics of torsional Alfven 
waves can be considered by adding non-linear 
terms to Eqs. (5)-(6), which then can be rewritten 
as 



<%0 

at 



dVs 



d 
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dZ dZ 
^[V R b 4 



V^b R ), 



(22) 



'dV+ VR d(RVt) 
{PQ + Pl) ^ + ll^R- 



V z 



dVt 
dZ 
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(B a + b z ) db^ b R d{Rb^) 
H dZ Rn dR 



(23) 



It is clear that the action has a parametric ori- 
gin, because all non-linear terms include at least 
one component of the torsional waves. Here we 
assume that the energy of slow magnetoacoustic 
waves is much larger than the energy of initial tor- 
sional waves. Then the back-reaction of torsional 
waves can be neglected at the initial stage and the 
slow magnetoacoustic waves can be considered as 
pumping waves with constant amplitudes. 

Recently IZaqarashvili fe Roberts! (2006) sug- 
gested that sound and linearly polarized Alfven 
waves are coupled when they propagate with the 
same speed along the magnetic field (i.e. f3 ~ 1) 
and the frequency (as well as the wave number) of 
the sound waves is twice the frequency of Alfven 
waves. In that case the action of the sound waves 
had a similar parametric origin. So a similar phe- 
nomenon can occur here, as in the high plasma 
limit the phase speeds of slow magnetoacoustic 
and torsional Alfven waves have similar values. 

In equations (5)- (6) the solution for torsional 
modes has an arbitrary i?-dependence. Therefore 
in the case of parallel propagating slow waves we 
may take V$ and to be independent on R in 
equations (22)-(23), which results in 

db^ dV<i, d d 

-dT = Ba ^z + dz {v * bz) -dz {Vzb * ) - 



(po + Pi) 



URb^ ^ hjiVtp 
Ro 
(dV 6 



\ dt 



Ro 
Ro 



V z 



dV* 
dZ 



(B + b z ) <%0 hub^ 

dz 



(24) 



(25) 



fi o/j Rq/j, 

The physics of wave coupling will be similar for 
obliquely propagating slow magnetoacoustic waves 
as well. In this case the influence of slow waves 
on the torsional waves can be considered at any 
cylindrical surface located at different R. Also the 
phase speed of slow waves along the magnetic field 
coincides with the Alfven speed, thus the resonant 
conditions given in IZaqarashvili fc Roberts! (|2006l ) 
will be satisfied. But, for simplicity, here we con- 
sider that both slow and torsional waves propagate 
along the magnetic field. 



In order to solve equations (24)-(25) we use the 
fo rmalism of weakly non-linear in teraction similar 
to Zaqarashvili fc Roberts! (|2006l ). So we consider 
the perturbations much smaller than the unper- 
turbed values. Then the variation of the torsional 
wave amplitudes due to the action of the slow 
waves will be a slow process. Therefore the per- 
turbed physical quantities of torsional waves can 
be represented as the product of a slowly varying 
amplitude and a rapidly oscillating term, 

= u^ty^'V + u%{t)e- l ^ z \ (26) 

& = k(t)^^ + %{t)e- ; ^ z \ (27) 

while the amplitude of slow waves remains con- 
stant 



(28) 



Vz 



u z e 



-iil> 2 (t,z) 



v>r = URe r2V ' } + u R e 



z< ■ (29) 



Here 



ipi = -Ld A t + k A Z, 
ijj 2 = -uit + k z Z, 



(31) 

(32) 
(33) 



where u>, k z , loa, are the frequencies and wave 
numbers of slow magnetoacoustic and torsional 
Alfven waves respectively. The frequencies and 
wave numbers satisfy the dispersion relation (8) 
which is the same for both waves. 

Substitution of expressions (26)-(31) into equa- 
tions (24)- (25), and averaging over rapid oscilla- 
tions in phase, leads to the disappearance of all 
exponential terms, so that only the first and third 
order terms remain. In the first approximation 
(neglecting all second and third order terms), the 
slow magnetoacoustic and torsional Alfven waves 
are decoupled and the amplitudes are constant. 
Third order terms (which are due to the advective 
terms in the momentum equation) are significant 
only in the case of very large amplitudes and so 
can be ignored. However, when 



V> 2 = ±2V>i, 



(34) 



the second order terms also remain after averaging 
and consequently the torsional wave amplitudes 
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become time-dependent. So if the frequencies and 
wave numbers satisfy the conditions 



lo = 2u A , k z = 2k A , 
u = -2u A , k z = -2k A , 



(35) 
(36) 



the waves are resonantly coupled and the slow 
magnetoacoustic waves may transfer energy to tor- 
sional Alfven waves. 

According to condition (35), averaging equa- 
tions (24)- (25) leads to the equations which gov- 
ern the temporal behaviour of the complex ampli- 
tudes: 



-iu A b^ + iuj A -^-b* <j} -ik A B {) u lt> = Q, (37) 

du<p f>i dil % ■ ~ . 2uij^ ik A B ot 
at po at R v jj,po 



2k A B 



u R b; = 0. 



(38) 



pp uj A R 

Here p\ = pn + ipu and ur = uri + iur2 are 
constants. 

Now let us search for solutions of (37)- (38) in 
the exponential form 



b<t> = (&0i + ib^e 81 , 



(39) 
(40) 



where b^i, b^, u^i and u^2 are constants. In ex- 
pressions (39)- (40) S is real, otherwise equations 
(37)-(38) are not satisfied. A positive S implies 
amplification of the initial perturbation, while a 
negative S implies damping. 

We substitute expressions (39)-(40) into equa- 
tions (37)-(38) and, without loss of generality, as- 
sume pn = u R2 = 0. Finally, we get equations 



5-u> A I &0i + cu A b^2 + kjiBoii^ = 0, (41) 

Po 



U A Uc/,1- 



k A B 0i 2k A B . ? . . 

-b(j>i 5-Ufli&02 = 0. (44) 



p,p uj A R 



(45) 



This system has a non-trivial solution if the sys- 
tem determinant is zero, which gives a fourth order 
equation for S. After dropping higher order terms, 
we get 

_ LO A Pl2 

2 po ' 

Thus the torsional Alfven waves grow exponen- 
tially in time and their growth rate depends on the 
relative amplitude of slow magnetoacoustic waves 
i.e. P12//O0, which is typical of weakly non-linear 
processes. 

In order to check the analytical solution, we 
performed the numerical simulation of spatially 
averaged equations (24)-(25) with k z = 2k A , so 
that only the time dependence is retained. Then 
the numerical calculation is straightforward and 
we find that torsional waves have an exponen- 
tially growing solution when their period is twice 
the period of the slow magnetoacoustic waves i.e. 
u> = 2ui A (Fig. 2). The up-plots refer to the density 
and velocity components of the slow waves (pi and 
ur) and the down-plots refer to the components 
of the torsional waves (b^ and u$). Here is as- 
sumed to be ~ 10 2 and the strong amplification of 
torsional waves begins after ~ 10-15 wave periods. 

Note that the growth rate of torsional Alfven 
waves depends on the value of plasma 0. The 
energy transfer occurs only when the density per- 
turbations in slow magnetoacoustic waves are 
nonzero. When goes to infinity, then the den- 
sity perturbations vanish and so there is no energy 
exchange between the waves. Thus the resonant 
process is of importance for large, but finite 0, 
which requires strong magnetic fields in the stel- 
lar interiors. 

5. Discussion 



S + u) A — ) - wa&^i - k A B u^i = 0, (42) 
Po J 

2k A B „ - k A B Q ~ 

-Mflittyi H = 0, (43) 



p,p a u} A R 



PPo 



We have seen that slow magnetoacoustic and 
torsional Alfven waves are resonantly coupled 
when they propagate with the same Alfven speed 
along an unperturbed magnetic field in a 3> 1 
medium. The slow magnetoacoustic wave para- 
metrically triggers the torsional wave with the half 
frequency and wave number. The growth rate of 
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the torsional mode depends on the amplitude of 
the slow magnetoacoustic wave: the larger is the 
amplitude of the slow magnetoacoustic wave, the 
faster is the energy transfer to torsional waves, 
which is a common scenario for non-linear pro- 
cesses. The growth rate depends on the value 
of plasma (3 as well; higher (3 leads to weaker 
coupling. The reason is that the density pertur- 
bations in slow magnetoacoustic waves are weaker 
for higher (3 and thus the slow waves do not alter 
the Alfven speed. However the situation probably 
will be changed for nonadiabatic slow waves as 
the density perturbations can be larger due to the 
compensation by the temperature variation. Then 
the density perturbations will affect significantly 
the Alfven speed, thus the torsional waves. This 
process seems to have interesting consequences 
and will be studied in future. 

The coupling between these waves may have 
interesting applications in astrophysics. The ex- 
istence of torsional oscillations in stellar inte- 
riors is a long time dilemma. Hydromagnctic 
torsional oscillations of a seed magnetic field 
have been suggested to take place in stellar in- 
teriors, which may lead to the cyclic ch ange in 
sign o f the sol a r toro i dal magnetic field dWalen 
1949[ ICowlind Il953t IPiddinetonl Il97lt iGough 
19881 ). However the oscillations probably are 
damped because of the transport of the mag- 
net icJhix_to_th£_sjarj^^ phase mix- 
ing ( Charbonneau fc MacGregorll993l ). Therefore 
the absence of a plausible energy source which 
can support torsional oscillations has been an ar- 
gume nt against their existence (jRosner fc Weiss 
1992). However, several different sources have 
previously b een suggested to support torsional 
oscillations (lLavzer et al Il979i : iMcIntvrd Il994 ; 
Lanza et al. I ll998t lGoughT l2002). but they are not 
convincing. 

Recently stellar radial pulsation in the funda- 
mental mode, which can be considered as a stand- 
ing fast magnetoacoustic wave in the presence of 
magnetic field, has been sug gested as a source 



for t orsional Alfven waves (jZaqarashvili et al 
2002). However, amplified torsional waves can 
not give rise to torsional oscillations due to 
their short wavelength as compared to the stel- 
lar radius (because the Alfven speed is much 
smaller than the sound speed). Therefore, rather 
than give rise to large-scale torsional oscillations, 



they will propagate upwards to the stellar atmo- 
sphere, and can be observed in the stellar wind 
( Zaqarashvili and Belvederdl2005l ). 

It is a natural physical process that in closed 
systems (the Sun or a simple tuning fork) almost 
all energy is stored in the fundamental mode. For 
example, take an ordinary tuning fork: any force 
(whatever its origin is) leads to an acoustic oscil- 
lation in the fundamental frequency. Obviously, 
there are other overtones, but with much smaller 
amplitudes. However, if one puts the tuning fork 
in a magnetic field embedding a plasma, then the 
acoustic frequency will split into two frequencies 
corresponding to fast and slow magnetoacoustic 
waves. Similarly, the large-scale magnetic field 
may cause the splitting of the fundamental acous- 
tic mode into high and low frequency branches, 
corresponding to fast and slow magnetoacoustic 
waves. Therefore a small deviation of the star 
from equilibrium must lead to the oscillation in 
both, fast and slow fundamental modes. Due to 
the small phase velocity, the period of standing 
slow magnetoacoustic waves will be much longer 
than the period of the fast fundamental mode. 
In the case of the adopted simple configuration 
of poloidal magnetic field in a stellar interior (see 
Fig.l), the standing slow magnetoacoustic waves 
can be pictured as radial oscillations with veloc- 
ity nodes along the Z axis, which causes a peri- 
odic variation in stellar oblateness The number of 
nodes depends on which harmonic is excited. As 
the phase speed is Alfvenic, the first harmonic will 
have a period of order 2Rq / ca, the second of order 
Rq/ca etc. Due to the radial polarization of the 
velocity, the oscillations may easily gain energy 
from various sources like gravity, pressure gradi- 
ent, rotation etc. 

Because of the parametric coupling, these oscil- 
lations may transfer energy to the torsional Alfven 
waves with twice the period and wavelength (see 
equation (35) and Fig. 2). The amplified torsional 
Alfven waves may easily lead to the set up of tor- 
sional oscillations as their wavelength now is com- 
parable to the stellar radius. Torsional oscillations 
now can be sustained against damping, as they 
may continuously gain energy from the radial pul- 
sation. 

The energy that may support the slow magne- 
toacoustic fundamental mode can be large. With 
reference to subsection 5.1., the gravitational en- 
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ergy of the Sun (i.e. the work required to dissi- 
pate solar matter to infinity), with radius Ro ~ 
7-10 10 cm and mass M w 2-10 33 g, is 6.6-10 48 
erg, even larger than the total in ternal radiant 
energy ps 2.8-10 47 erg (jAllenl Il992h . The global 
density variation due to the slow magnetoacoustic 
fundamental mode is accompanied by a gravita- 
tio nal potential variation ( see the classical paper 
bv lLedoux fe Pekerid ( 194~ih ). which unfortunately 
is neglected in this paper due to the simplicity 
of presentation. Suppose that the radius varia- 
tion AR due to the slow magnetoacoustic fun- 
damental mode is as small as AR/Rq ~ 10~ 5 
(for the Sun it is ~ 0.01", however observations 
also show radius variations ~ 0.02" dEmilio et al.1 



l2000h and even - 0.08" (jLaclare et all Il996)). 



The gravitati onal energy change for a homoge- 
neous sphere ( Chandrasekhai 19391 ) will be of the 
order of 6.6-10 48 Ai?/i? ~ 6.6-10 43 erg (this is the 
energy stored in the fundamental mode). Sup- 
pose that only 0.1% of this energy goes into the 
energy of torsional oscillations. Then the cor- 
responding magnetic energy variation will be of 
the order of ~ 6.6T0 40 erg. Even if the mag- 
netic field is redistributed throughout the whole 
stellar volume, we get a toroidal magnetic field 
variation with an amplitude > 10 4 G. So the en- 
ergy is more than enough for supporting the tor- 
sional o scillation against ohmic diffusion o r phase 
mixing ( Charbonneau fc MacGregorlll993 ). More- 
over, the amplification of the toroidal magnetic 
field component will be accompanied by the mag- 
netic buoyancy instability, which leads to the erup- 
tion of magnetic flux towards the surface. Thus 
torsional oscillations can not cause a significant 
back-reaction on the pumping radial pulsation as 
their further amplification will be balanced by the 
energy transport towards the surface in form of 
magnetic flux tubes, which leads to chromosphcric 
and coronal activity. 

However, it must be mentioned that the many 
simplifications of the model (e.g. uniform density, 
cylindrical geometry) will cause significant limita- 
tions to the application of the mechanism to real 
situations in stars. For instance, notice that the 
density (and the magnetic field) can be considered 
as uniform only in a narrow cylindrical shell, not 
throughout the star, but the physics of the mech- 
anism still remains substantially the same. Con- 
sidering spherical symmetry instead of cylindri- 



cal one will be a further improvement, however it 
will cause additional mathematical complications. 
The aim of this paper is to show the existence of 
such an energy transformation mechanism, which 
can be developed and applied to more realistic as- 
trophysical situations in future. 

Here we briefly discuss possible astrophysical 
applications of the wave coupling discussed in this 
paper. 

5.1. Solar radius variation and torsional 
oscillations 

The variation of the solar radius during the 
11-year period activity cycle is a long time 
problem. The characteristics of the variation 
may reveal some basic mechanisms of activity, 
therefore many efforts have been made in order 
to find observational evidence of a correlation 
between radius variation and activity parame- 
ters. The observational results are controver- 
sial and somewhat inconsistent. Some observa- 



tions indicate a negative correlation (jLaclare et al 



19961 ) . so me results show no systematic variatio n 
in time ( Brown fc Christensen-Dalsgaard 1998). 
Most observations reveal a positive correlation 
( Ulrich fc Bertellolll995t iBasul Il998l : lEmilio et al 



2000; iNoell l2004h . The observed amplitude of 
the variation is also controversial. Ba s ed on 
Mount Wilson data llllrich fc Bertellol (| 19951 ) 



found ~ 0.2"; INoell (|2004l ) suggested a s imilar 
amplitude ~ 0.2", while in lLaclare et al. (1996) 
the amplitude is ~ 0.08". Probably, the most real- 
isti c value of the ampl itude ~ 0.02" was obtained 
by lEmilio et all (|2000h by analyzing SOHO-MDI 
data. The variation of the solar radius during the 
activity cycle naturally leads to the idea of its con- 
nection to the magnetic field pressure variation in 
the interior due to the large-scale turbulent dy- 
namo. However, even a very small variation of the 
radius requires a huge magnetic energy to com- 
pensate the corresponding gravitational energy 
change (see the above estimate). Then it seems 
plausible that the radius variation is due to the 
global slow magnetoacoustic mode. Slow magne- 
toacoustic waves propagating along the poloidal 
magnetic field have both R and Z-components 
of velocity (see Eqs.(9)-(14)), therefore the global 
slow mode will lead to the periodical variation of 
the solar oblateness, which may cause interesting 
consequences for solar total irradiance. The fun- 



S 



damental period of pulsations in the global slow 
magnetoacoustic mode will be of order 2Rq/ca, 
which for the 11-year period implies the mean 
magnetic field strength in the solar interior to be 
~ 10 3 G. 

Another possible candidate for the global mode 
coupling can be the observed 1.3 yr variation at 



the b ase of the solar convection zone (jHowe et al 
2000) . If the global slow mode has the period of ~ 
0.65 yr, then it may drive the torsional oscillations 
with the period of 1.3 yr in the solar tachoclinc. 

5.2. Orbital period modulation in close bi- 
naries 

Long-term observations of close binary stars 
show a variation of the orbital period over long 
time-scales with relative amplitude AP/P ~ 
10~ 5 , where P is the orbital period. A variety 
of theoretical models have been proposed to ex- 
plain this phenomenon. The first models included 
the observed period changes due to the apsidal 
motion, or the i nfluence of a dista nt, unseen com- 
panion. Then Applegatel ( 19921 ) suggested that 
the variation can be explained by the gravitational 
coupling of the orbit to the variation in shape of 
a magnetically active companion. The variable 
deformation of the active star is produced by the 
variation in the distribution of angular momen- 
tum during t he ac tivity cycle. On another hand, 
Lanza et al. (1998) suggested that a torsional os- 
cillation driven by the stellar dynamo may account 
for the observed behaviour of the binary system. 
However the problem still remains, because the 
magnetic energy can not compensate the gravita- 
tional energy change, as has been shown for the 
Sun. 

Instead, here we propose that the long term ra- 
dial pulsation due to the global slow mode consid- 
ered in our paper (which causes a periodic vari- 
ation in oblateness) may be responsible for the 
variation of the orbital period with a gravitational 
coupling similar to Applegate's. However contrary 
to Applegate's mechanism, here stellar radial pul- 
sation can be the reason for both orbital period 
modulation and magnetic activity in the primary 
star. 

The relative amplitude of the orbital period 
modulation is proportional to the relative varia- 
tion of the gravitational acceleration of the pri- 



mary AP/P ~ Ag/g (|Applegatejll992f ). On the 
other hand, the relative variation of the gravi- 
tational acceleration is of the order of Ag/g ~ 
AR/R, where R is radius of the primary. Then in 
order to get the observed relative period modula- 
tion AP/P ~ 10 -5 , the radius of primary needs 
to vary as much as AR/R ~ 10~ 5 . It is very 
intriguing that the observed solar r adius relative 



variation deduced from SOHO data ( Emilio et al 



2000) has a comparable extent AR/R ~ 2-10" 5 . 
So relative radius variations of the primary with 
amplitude similar to the Sun's, can explain the 
observed orbital period modulation in binary sys- 
tems. 

6. Conclusion 

In this paper we considered the weakly non- 
linear coupling between slow magnetoacoustic and 
torsional Alfven waves in the high plasma (3 limit. 
It has been shown that slow magnetoacoustic 
waves propagating along the unperturbed mag- 
netic field may resonantly amplify torsional Alfven 
waves with the half frequency and wave-number of 
slow magnetoacoustic waves. The growth rate of 
torsional Alfven waves depends on the amplitude 
of slow magnetoacoustic waves and the value of 
plasma j3. Therefore the global slow magnetoa- 
coustic waves, excited due to the action of various 
forces (gravity, gas and radiation pressure gra- 
dients, rotation etc.) in a stellar interior, may 
give rise to large-scale torsional oscillations. The 
process may be of importance in solar and stellar 
physics, but needs more detailed study. 
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Fig. 2. — The time evolution of spatially averaged 
(with Z) slow magnetoacoustic (up) and torsional 
modes (down). The strong amplification of the 
torsional mode with twice the period of the slow 
magnetoacoustic mode is clearly seen, as predicted 
by the analytical solution. Here p%/po — 0.001 is 
considered. 



12 



